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1 Introduction 



Atmospheric muons observed in the deep atmosphere are carrying informa- 
tion on the interaction of primary and secondary cosmic-ray particles in the 
atmosphere. Since decay probabiUty and energy loss rate depend on energy, 
high-statistics measurement of spectral shape can provide useful knowledge 
to understand propagation process of cosmic-ray particles through the atmo- 
sphere. For a more detailed study of neutrino oscillation phenomena observed 
in atmospheric neutrinos it is essentially important to reduce the systematic 
errors in the predicted energy spectra of neutrinos. In order to improve the 
accuracy of the predictions, a detailed understanding of hadronic interactions, 
as well as precise information of primary cosmic-ray fluxes, is indispensable. 
As for absolute fluxes of cosmic-ray protons and hehum nuclei below 100 GeV, 
a very precise measurement was carried out by using the BESS detector [1]. 
Two independent experiments, the AMS [2,3] and the BESS, show extremely 
good agreement with each others in their results. Since production and decay 
process of muons are accompanied by neutrino productions, the muons fluxes 
correlate directly to the hadronic interactions. Thus the precise measurement 
of atmospheric muon spectra can also provide important information to im- 
prove accuracy of atmospheric neutrino calculations. 

We measured absolute fluxes of atmospheric muons at the top of Mt. Norikura, 
Japan with the BESS detector, which was the same apparatus as we utilized to 
measure sea-level muons [4] as well as the primary protons and helium nuclei 
[1]. These measurements provide fundamental and very useful information 
about hadronic interactions and geomagnetic effects. 



2 BESS experiment 

2.1 BESS Detector 

The BESS (Balloon-borne Experiment with a Superconducting Spectrometer) 
detector [5,6,7,8,9] is a high-resolution spectrometer with a large acceptance to 
perform highly sensitive searches for rare cosmic-ray components, and precise 
measurement of the absolute fluxes of various cosmic rays. Fig. 1 shows a 
schematic cross-sectional view of the BESS instrument. In the central region, a 
uniform magnetic held of 1 Tesla is provided by using a thin superconducting 
solenoidal coil. A magnetic-rigidity {R = Pc/Ze) of an incoming charged 
particle is measured by a tracking system, which consists of a jet-type drift 
chamber (JET) and two inner- drift- chambers (IDC's) inside the magnetic fleld. 
The deflection (i?^^) and its error are calculated for each event by applying 
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a circular fitting using up-to 28 hit points, each with a spatial resolution 
of 200 fim. The maximum detectable rigidity (MDR) was estimated to be 
200 GV. Time-of-flight (TOF) hodoscopes provide the velocity and energy 
loss {dE/ dx) measurements. A 1/P resolution of 1.6 % was achieved in this 
experiment. 

Particle identification was made by requiring proper 1/(3 as well as dE/ dx 
as a function of the rigidity. An electromagnetic shower counter has been 
equipped for electron/muon separation. It consists of a 2 thick lead plate 
and an acrylic Cerenkov counter. The lead plate covers about 1/5 of the total 
acceptance. 

The simple cyhndrical shape and the uniform magnetic field make it simple 
and reliable to determine the geometrical acceptance precisely. The hve data- 
taking time was measured exactly by counting 1 MHz clock pulses with a 
scaler system gated by a "ready" status that controls the first-level trigger. 
The resultant live-time ratio was as high as 98.8 % through the measurement. 

2.2 Data samples 

The atmospheric cosmic-ray events were observed at Norikura Observatory, 
ICRR, the University of Tokyo, Japan. It is located at 36° 06' A^, 137° 33'^;, 
which is about 40 km southeast of Kamioka. The altitude is 2,770 m above sea 
level. The vertical geomagnetic cutoff rigidity is 11.2 GV [10]. The experiment 
was performed during two periods of 17th - 19th and 21st - 23th of Septem- 
ber 1999. During the observation, the atmospheric pressure and temperature 
varied as shown in Fig. 2. The mean (root- mean-square) atmospheric pressure 
and temperature were 742.4 (2.9) g/cm^ and 10.9 (1.1) °C, respectively. The 
total live data-taking time was about 4 days. 



3 Data analysis 

3.1 Data Reduction 

In this analysis, the BESS detector was divided into two sections as shown 
in Fig. 1; i.e., the "normal" section and the "EM shower" section. The EM 

shower section includes the shower counter. The normal section was used to 
measure muon momentum spectra. In order to estimate small contaminations 
of electrons, the data samples recorded in the EM shower section were exam- 
ined. 
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In the first stage of data reduction, we selected events with a single-track 
fully contained within a fiducial volume. The fiducial volume was defined as a 
central region in the JET chamber. This definition ensured the longest track 
to achieve the highest resolution in the rigidity measurement. The zenith an- 
gle {9z) was limited to be cos^^; > 0.98 to obtain near vertical fluxes. The 
single-track event was defined as an event which has only one isolated track 
and one or two hit-counters in each layer of the TOF hodoscopes. Since the 
major component of the observed particles was muons, those which can easily 
penetrate the detector, almost all particles passed this selection criterion. The 
single-track selection eliminated very rare interacting events. In order to ver- 
ify the selection, events were scanned randomly. It was confirmed that 1,995 
out of 2,000 visually-identified single-track events passed this selection criteria 
and interacting events were fully eliminated. Thus, the track reconstruction 
efficiency was determined to be 99.8 ± 0.1 % for the single-track event. The 
single-track events selected in this stage were followed by the succeeding anal- 
ysis. 



The dE/ dx inside the TOF hodoscope was examined to ensure that the parti- 
cle was singly charged. The particle identification was performed by requiring 
proper l/f3 as a function of the momentum as shown in Figs. 3 and 4. Since 
the 1//9 distribution was well described by Gaussian and a half width of the 
1/(3 selection band was set at 3.89 a, the efficiency was very close to unity. 
Taking into account a small deviation from pure Gaussian, the efficiency was 
evaluated to be 99.9 %. In order to assure accuracy of the momentum measure- 
ment, event qualities such as in the track fitting procedure were required. 
The efficiency of this quality-check was evaluated to be 96.3 ±0.1 % as a ratio 
of the number of good quality events to all events between 10 and 20 GeV/c. 
Since the quality-check required consistency between the fitted track and each 
hit in the tracking detector, its efficiency does not depend on the momentum 
in higher momentum region where multiple scattering can be neglected. Ac- 
cordingly the above value of the quality-check efficiency was applied to higher 
momentum region. Before and after the quality-check and succeeding correc- 
tion of its efficiency, the obtained flux above 10 GeV/c changed no more than 
±2.2 %. These changes were treated as systematic errors in the quality-check. 
Below 10 GeV/c, quality-check was not imposed, because momentum reso- 
lution before the quality-check was high enough to measure the momentum. 
Thus there was no systematic error associated with the quality-check. 



With the data reduction described above, 707,563 muon candidates were se- 
lected. The combined efficiency in the data reduction process was 99.6 ±0.1 % 
and 95.8 ± 2.2 % below and above 10 GeV/c, respectively. 
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3.2 Contamination estimation 



As seen in Figs. 3 and 4, main background sources were protons for positively 
charged muons whereas electrons and positrons for both positively and nega- 
tively charged muons. Protons were identified by examining l/f3 distribution. 
Below 1.6 GeV/c, protons were clearly separated from muons. Between 1.6 
and 3.2 GeV/c, the proton contamination in the muon candidates and its 
error were estimated by fitting the 1/(3 distribution with a double-Gaussian 
function as shown in Fig. 5. The ratio of observed protons to muons was found 
to be well described as a simple power law as p/fi~^ = 0.24 • p-o-^o in a mo- 
mentum region of 1 - 3 GeV/c. The proton contamination above 3 GeV/c was 
estimated by extrapolating this power law. Above around 10 GeV/c, however, 
the pZ/x"*" ratio will be constant. It is reasonably assumed that pions are pro- 
duced having a similar spectral index of the parent protons and then decay to 
muons, those which can reach the mountain altitude before they decay in a 
higher momentum region. Under this assumption, the spectral index of muons 
is similar to that of the parent protons, resulting in a constant p/yU"*" ratio. 
This tendency was well reproduced by two Monte Carlo simulations [11]; one 
simulation incorporated the FRITI0F7 [12] as a hadronic interaction model 
and the other was developed based on the DPMJET-III event generator [13]. In 
this analysis, the p/yU"*" ratio above 10 GeV/c was assumed to take a constant 
value of 3.8%. These proton contaminations were subtracted from the muon 
candidates. An error originated from the proton subtraction between 1.6 and 
3.2 GeV/c was estimated during the fitting procedure to be less than 1 %. 
In a higher momentum region, the p/n'^ ratios predicted by the two Monte 
Carlo simulations agreed with each other within ±1%. These discrepancies 
were included into systematic errors. 

The muon candidates might include pion contamination, since it has the same 
signature and is hardly identified. According to the same Monte Carlo simula- 
tions that was utilized for p/ fi^ ratio estimations, tt/ fi ratio was estimated to 
be less than 1 % below 40 GcV/c and less than 3 % between 40 and 100 GeV/c. 
It is smaller than the statistical errors over the whole momentum range in this 
analysis. The pions were neglected in this analysis both in contamination sub- 
traction and systematic error estimation. 

The electron contamination was to be considered as a background well below 
1 GeV/c. For an estimation of electron contamination, events passed through 
the EM shower section were considered. The Cerenkov light yield from the 
shower counter was examined to estimate an e//x ratio. Fig. 6 demonstrates 
the following estimation procedure. The Cerenkov light yields for electrons and 
muons were individually obtained from Monte Carlo simulations. The observed 
light output distribution is a sum of these two distributions with appropriate 
weights. By changing weights for those two distributions, the best fit e/// ratio 
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was obtained. The electron contamination in the muon candidates observed in 
the EM shower section was about 1.5 % at 0.6 GeV/c and rapidly decreased 
with increasing momentum. The e/fi ratio below 1 GeV/c was found to be 
described as e/// = 0.0024 • P~^-^. Above 1 GeV/c, the number of observed 
electron events was too small to estimate its ratio by this procedure. We esti- 
mated the electron contamination by extrapolating the ratio obtained below 
1 GeV/c to a higher momentum region. After correcting the different data 
reduction efficiency in the normal section from that in the EM shower section 
by using the Monte Carlo simulations, the electron contamination was sub- 
tracted from the muon candidates observed in the normal section. Accuracy 
of the electron subtraction was limited by poor statistics of electron events. 
An error in this subtraction was estimated from the fitting procedure to be 
1.0 % at 0.6 GeV/c and rapidly decrease with increasing momentum. 

The proton contamination was to be properly subtracted above a few GeV/c. 
On the other hand, the electron contamination should be carefully considered 
and properly subtracted in a different momentum region below 1 GeV/c. Thus 
the overall error in the contamination subtraction was as small as 1 % over 
the whole momentum range. It may be worth mentioning that e//i and p/ n 
ratios discussed above are different from flux ratios because the efficiencies for 
protons and electrons are not same as that for muons. 



3.3 Corrections 

In order to determine the atmospheric muon spectra at the mountain altitude, 
the following corrections are required; (i) exposure factor, (ii) ionization energy 
loss, and (iii) interaction loss. 

The exposure factor is a product of geometrical acceptance and live-time. The 
geometrical acceptance defined for this analysis was calculated by simulation 
technique [14] to be 0.0214 m^sr for energetic particles which have straight 
track inside the tracking volume. It was almost constant over the whole mo- 
mentum range. The simple cylindrical shape and the uniform magnetic field 
make it simple and reliable to determine the geometrical acceptance precisely. 
The error arose from uncertainty of the detector alignment was estimated to 
be 0.4 %. Between observed and simulated zenith angular distributions of the 
trajectories, a small discrepancy was found below 2 GeV/c. In the simulation 
isotropic trajectories were generated to obtain the geometrical acceptance, 
however, the flux of atmospheric muons has a zenith angle dependence. In 
a low momentum region, the dependence might be significant even in the so 
near vertical direction as cos ^2 > 0.98. The discrepancy was included in the 
systematic errors. The systematic error associated with this discrepancy was 
estimated to be 3.5 % at 0.6 GeV/c and decreased with increasing momentum. 
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The live data-taking time was measured exactly to be 339,874 seconds. 

The energy of each particle at the top of the instrument was calculated by 
summing up the ionization energy losses inside the instrument with tracing 
back the event trajectory 

In order to estimate the interaction loss probabihties inside the BESS detec- 
tor, Monte Carlo simulations with the GEANT code [15] were performed. The 
probabilities that muons can be identified as single-track events were evalu- 
ated by applying the same selection criterion to the Monte Carlo events. The 
resultant single-track efficiency was 97.9 % at 10 GeV/c and almost constant 
over the momentum region discussed here. The Monte Carlo simulation of the 
BESS detector well reproduced the observed event profile. By comparing the 
observed and simulated parameters used in the single-track selection, i.e., the 
number of tracks inside the tracking volume and the number of hit-counters 
in each layer of the TOF hodoscopes, it was concluded that a discrepancy 
between the observed and simulated event shapes was no more than ±1.6% in 
the whole momentum range discussed here. This discrepancy was included in 
the systematic error of the single-track efficiency. Another source of the sys- 
tematic error in the Monte Carlo simulation was imccrtainty of the material 
distribution inside the BESS detector, which was estimated to be ±10 %. Since 
the single-track efficiency was as high as 98 %, the systematic error originated 
from this uncertainty was as small as 0.2 %. 

The combined systematic errors obtained by quadratically adding the uncer- 
tainties were ±3.0, ±1.4, and ±1.7 % at 1, 10, and 100 GeV/c, respectively. 



3.4 Spectrum deformation effect 



Because of the finite resolution in the rigidity measurement and the steep 
spectral shape, the observed spectrum may suffer deformation. In a low rigidity 
region well below the MDR, the spectrum deformation is negligibly small. In a 
higher rigidity region, however, a deflection becomes comparable to an error in 
curvature measurement; thus, the spectrum deformation is to be considered. 
This spectrum deformation effect of the BESS detector was discussed in our 
previous paper [1]. According to the discussion, the effect was estimated to be 
no more than 3 % below 120 GeV/c. The deformation effect is smaller than 
the statistical errors over the whole momentum range in this analysis. Thus 
no correction was applied to the measured spectra. 
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4 Results and discussions 

The momentum spectra of the atmospheric muons in momentum ranges of 0.6 
through 106 GeV/c at a mountain altitude of 2,770 m above sea level have 
been measured with the BESS detector. The results are summarized in Table 1 
and shown in Fig. 7. The first and second errors in Table 1 represent statistical 
and systematic errors, respectively. The overall errors including both errors are 
less than ±10 %. In Fig. 7, our result is compared with other absolute flux 
measurements [4,16,17]. A discrepancy between the observed muon fluxes is 
mainly due to the difference in altitudes. 

Fig. 8 shows charge ratios of atmospheric muons observed in a series of muon 
measurements with the BESS detector. The big difference in the charge ratios 
observed in Japan and Canada comes from the influence of the geomagnetic 
cutoff rigidity as discussed in Ref. [4]. Below 10 GcV/c, the charge ratio at 
mountain altitude is slightly lower than that at sea level. This tendency was 
made by energy loss inside the atmosphere between mountain altitude and sea- 
level. In a higher momentum region, the ratio seems to approach to around 
1.28. 

Our results were compared with theoretical predictions, in which atmospheric 
muon spectra were calculated simultaneously with atmospheric neutrinos. 
Fig. 9 shows comparisons between our results and the predictions. The pre- 
dicted spectra were calculated based on Ref. [18], but some parameters, such 
as atmospheric pressure, temperature and zenith angle limit, were adjusted 
to reproduce the experimental conditions. In Ref. [18], two calculations were 
discussed; one was developed on the basis of the HKKM interaction model 
[19] and the other was calculated by employing the DPMJET-III event gener- 
ator [13] as hadronic interaction models. As Fig. 9 indicates, the DPMJET-III 
interaction model reproduced the observed spectra better than the HKKM 
interaction model. Some deviation, however, was found between the observed 
and calculated spectra especially below 1 GeV/c. It seems reasonable to sup- 
pose that the hadronic interaction models are to be tested and improved by 
accurate spectra of atmospheric muons. The precise measurement of atmo- 
spheric secondary particles will improve the accuracy in the atmospheric neu- 
trino calculations. 



5 Summciry 

We have measured absolute fluxes of atmospheric muons at the top of Mt. 
Norikura, Japan. It is located at an altitude of 2,770 m above sea level. The 
overall errors of our measurements are less than ±10 %. The absolute fluxes 
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showed relatively good agreement with theoretical predictions calculated by 
using the DPMJET-III hadronic interaction model. The precise measurement 
of the atmospheric muons can modify the hadronic interaction model and 
improve accuracy of atmospheric neutrino calculations. 
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Fig. 1. Schematic cross-sectional view of the BESS instrument. 



11 



750 




17 18 19 20 21 22 23 

DATE (September 1999) 



Fig. 2. Atmospheric pressure and temperature during the observation. The experi- 
ment was performed during two periods of 17th - 19th and 21st - 23th. 
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Fig. 3. Positive muon selection in vs. momentum. Muons were selected with 
high efficiency. Proton and positron contaminations were subtracted. 
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Fig. 4. Negative muon selection in vs. momentum. Muons were selected with 
high efficiency. Electron contamination was subtracted. 
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Fig. 5. An estimation of proton contamination by TOF information. Muons have a 
peak at A(l//3) = 0. Filled circles are data at around 2.7 GeV/c. Best fitting curves 
are also shown. 
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Fig. 6. An estimation of electron contamination in events of the shower counter. 
Filled circles are data between 0.3 and 0.6 GeV/c. Histograms show the best fit 
Monte Carlo data. 
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Fig. 7. Absolute differential momentum spectra of atmospheric muons. Vertical axis 
shows the sum of and fi~ fluxes multiplied by P^. Only statistical errors are 
included. The geomagnetic cutoff rigidity at each site is shown as Re- The difference 
among the fluxes is mainly due to the difference in altitudes. 
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Fig. 8. Charge ratios of atmospheric muons. Only statistical errors are included. 
The geomagnetic cutoff rigidity at each site is shown as Re- The difference among 
the ratios is caused by the difference in geomagnetic cutoff rigidity and energy loss 
inside the atmosphere. 
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Fig. 9. Comparison between the observed and calculated spectra of atmospheric 
muons. The calculation with the DPMJET-III event generator reproduced our re- 
sults better than the HKKM interaction model. 
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Table 1 

Absolute differential momentum spectra of atmospheric muons. 
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